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' Abstract. Within the framework of chiral perturbation theory with virtual photons and leptons, we present 

an updated analysis of the pionic beta decay including all electromagnetic contributions of order e^'f? . 
We discuss the extraction of the Cabibbo-Kobayashi-Maskawa matrix element from experimental 

data. The method employed here is consistent with the analogous treatment of the Kiz decays and the 
determination of IKial- 
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1 Introduction 

The determination of the quark mixing matrix Vckm is 
presently a major issue in elementary particle physics. The 
existence of a direct CP violating phase in the kaon sec- 
tor has now been clearly established experimentally 
Furthermore, the failure of the CKM matrix elements to 
satisfy the constraints expressing the unitarity of Vckm 
would with certainty establish the existence of new de- 
grees of freedom beyond those present in the three gen- 
eration standard model. Since \yuh\ ^ 5 x 10~^ (the lat- 
est value obtained by CLEO from inclusive semileptonic 
B decays reads |Kf,| = (4.08 ± 0.63) x lO^^ g) and 
given the present accuracies on |Vud| and at the 

0.1% and 1% level, respectively, the most stringent test 
of the unitarity of the CKM matrix comes from the light 
quark sector, where + ^Yus^ should show compati- 

bility with unity with an accuracy better than 0.3%. The 
semileptonic K^z decay modes are presently considered to 
provide the best determination of with the current 
value \Vus{'K^z)\ = 0.2196 ± 0.0026 g, while the most 
accurate measurement of relies on the Tt values ob- 
tained from the super-allowed Fermi transitions of sev- 
eral 0+ nuclei, \Vud{J"C)\ = 0.9740 ± 0.0005 g]. This gives 
|Kd(^t)|2 + \Vus{K,^)? = 0.9969(17), i.e. a deviation 
from unity by 2a. It is of course too early to draw definite 
conclusions from this result. On the one hand, the value 
of \Vud\ obtained from the nuclear beta decays hinges on 
the control of the nuclear structure aspects involved in 
the evaluation of the radiative corrections to the transi- 
tion matrix elements. On the other hand, the extraction 



of \Vus\ from the Kf.3, decay modes relies on a one- loop 
chiral perturbation calculation [^||] of the relevant form 
factor, supplied with a model dependent estimate of the 
higher order contributions. The Particle Data Group com- 
pilation therefore recommends to double the uncertainty 
in the value of | Vud \ inferred from the nuclear Fermi tran- 
sitions, |Kd(.?^t)| = 0.9740±0.0010. In the case of the 
decay mode, other quark model estimates of higher order 
corrections lead to a wider range of values than those con- 
sidered in A similar conclusion, based on the analysis 
of 0{m1) tree- level corrections in both Fk/Ft^ and ifes 
form factors, has also been reached in [^[]. Determinations 
of I Kid I from other sources than the nuclear Fermi tran- 
sitions have therefore been considered. The neutron beta 
decay, where radiative corrections can be evaluated in a 
more reliable way, suffers from the drawback that in this 
case both the vector current and the axial-vector current 
contribute. It is thus not enough to measure the lifetime 
only, but the asymmetry in the electron emission angle 
with respect to the neutron polarization is also required. 
A recent measurement of this asymmetry [|| leads to a 
value |Kid('T-/3)| = 0.9713 ± 0.0013, even smaller than the 
one coming from the superallowed nuclear 0+ — > 0+ tran- 
sitions, and corresponds to \Vud{'n(i)\'^ + \Vus{Kes)\^ = 
0.9916(28), i.e. a ia deviation from unitarity. 

Another interesting possibility, which shares the ad- 
vantages of both Fermi transitions (pure vector transition, 
no axial-vector contribution) and neutron beta decay (no 
nuclear structure dependent radiative corrections) is pro- 
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expressions of the form factors have been 
but no evaluation of the relevant 0{p^) coun- 
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vided by the beta decay of the charged pion. The difficulty 
here Ues in the extremely small branching ratio, ^ 10~®. 
Nevertheless, such a measurement is presently being per- 
formed at PSI by the PIBETA collaboration , with the 
aim of measuring the branching ratio with 0.5% accuracy. 
At this level of precision, radiative corrections have to 
be taken into account, and the present paper addresses 
this problem. We shall use the effective theory formalism 
for processes involving light pseudoscalar mesons, photons 
and leptons introduced in pO| which is particularly well 
suited for the pionic beta decay. After a brief review of 
the main kinematic features of the process (Section 2), we 
describe the modifications to the structure of the decay 
amplitude induced by the radiative corrections in Section 
3. In particular, we obtain the corrections of order 0{e^p^) 
to the form factor. Section 4 discusses the radiative beta 
decay rate, which has to be included in order to cancel the 
infrared divergences which appear in the radiative correc- 
tions to the beta decay amplitude without real photon 
emission. Numerical estimates, in particular of the theo- 
retical uncertainty in the determination of |T4id(7''/3)|, are 
presented in Section 5. Conclusions are given in Section 
6, and two Appendices contain some technical material 
related to the calculation of the loop contributions. 



2 Kinematics 

In the absence of radiative corrections, the invariant am- 
plitude of the decay 



'^^{P+) 7i'"(Po) e+(pe) lyeiPi.) 



(2.1) 



reads 



-Po) 



where 



F = u{p,)Y{l~l5)v{pe)■ 



Po)^. 
(2.2) 

(2.3) 

The expression in parentheses corresponds to the matrix- 
element (7r°(po)|u7;id|7r"'"(p+))/-\/2. The hadronic form fac- 
tors depend on the single variable t = (p+ — po)^. 

The spin-averaged decay distribution p{y, z) is a func- 
tion of the two variables 



2p-(- ■ Po _ 



y = 



2p+ • Pe 



2E, 



(2.4) 



where E^^o {E^+ ) is the 7r° (positron) energy in the rest 
frame of the charged pion. Alternatively, one may also use 
two of the Lorentz invariants 



pof = A4+il + r^-z), 



where 



t = {p 

u={p+ -Pef = Ml+{l + r,-y) 
s^{po+Pef ^Ml+{y + z-l), 

ml _ Ml, 



(2.6) 



Then the Dalitz plot density (without radiative correc- 
tions) reads 



p^") {y,z)^N [Af^ I it) I ^ + /(°) (i)/l") it) 



with 



(2.7) 



AT 



(2.8) 



647r3 

r{Tr+ ^ Tr"e+iye) = j dydz (y, z). (2.9) 

V 

The kinematical densities are given by 

Af\y, z) = 4(z + y - 1)(1 - y) + r^i^y + 3z - 3) 
- 4r^ -I- re(r^ - re), 

4°^(y,z) = 2re(3-2y-z + re-r,), 



Af\y, z) = re(l + - z - re). 

The physical domain T) is defined by 

2\A^ < y < 1 + r-e - r^, 
a(j/) - b{y) < z < a{y) + b{y), 



(2.10) 



(2.11) 



where 

a{y) = 

Hy)-- 

or, equivalently. 



(2 - y) (1 + re + ?v - y) 
2(l + re-y) 



\Jy'^ - 4re (1 + re - r^ - y) 
2(1 + re -2/) 



(2.12) 



c(z) - d{z) <y < c{z) + d{z), 



(2.13) 



where 



c(z) = 
diz) = 



(2 -z){l + r^+re- z) 

2(l + r^-z) 
^/ z^ - 4r^ (1 + r^ - re - z) 
2(1 + r^ -z) ■ 



(2.14) 



3 Pionic beta decay and electromagnetism 

The presence of the electromagnetic interaction d oes not 
change the structure of the invariant amplitude (2.2) in 
terms of the form factors, but changes the form factors 
themselves [||: 

(2.5) M[/°\t)J^^\t)]^M[F+{t,u),F^{t,u)]. (3.1) 



The full form factors F± {t, u) contain the effects of vir- 
tual photon exchange and the contributions of the ap- 
propriate electromagnetic counterterms. These quantities 
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depend also on a second kinematical variable as they can- 
not be interpreted anymore as matrix elements of a quark 
current between hadronic states. 

Photon loop diagrams contributing to the weak vertex 
function are shown in Figure |^. In particular, it is the 
second diagram which introduces the dependence on the 
second kinematical variable u. 




7r+(p+) 



7r+(p+) 




7r+(p+) 



vrO(po) 



Fig. 1. Photon loop diagrams (wave function renormalization 
diagrams are not shown). The black dots denote (scalar) QED 
vertices, the box denotes the vertex proportional to Gf- 



The following observation considerably simplifies the 
further analysis of the pionic beta decay: the form factor 
F- (t, u) is proportional to the isospin suppressed mass dif- 
ference M^_|_ — M^o = 2e'^ZF^. In addition, the kinemat- 
ical densities A2 and ^3 m ultip lying in the formula 
for the Dalitz plot density ( ^.7|) are proportional to the 
small quantity ~ 10^^. Therefore, these contributions 
can safely be neglected and we may restrict ourselves to 
the discussion of F+. 

The form factor F-^-{t,u) contains infrared singulari- 
ties due to low-momentum virtual photons. They can be 
regularized by introducing a small photon mass M-y. The 
dependence on an infrared cutoff reflects the fact that 
F^{t,u) cannot be interpreted as an observable quantity 



but has to be combined with appropriate contributions 
from real photon emission to arrive at an infrared-finite 
result. 

It is convenient to decompose u) into a structure- 
dependent effective form factor f+{t), and a remaining 
part containing in particular the universal long distance 
corrections. Confining ourselves to electromagnetic contri- 
butions of order a, the full form factor is given by 

F+it,u)=[l + ^r{u,mlM^^;M^)\ f+{t). (3.2) 

Expressed in terms of the functions Fc, Fi, F2 defined in 
[D, F can be written as 

F{u, ml,Ml+-M^) = Fc{u, ml, M^^;M^) 
+ F,iu,mlM^+) 
+ F2{u,mlMl+). (3.3) 

The explicit expressions of the functions F^ A, -Hj can 
be found in Appendix]^. Fc{u, ml, M^) corresponds 
to the long distance component of the loop amplitudes 
which generates infrared and Coulomb singularities. (In 
our case, the Coulomb singularity lies outside of the phys- 
ical region.) The term Fi{u,ml,M^+) + F2{u,ml,M^+) 
represents the remaining nonlocal photon loop contribu- 
tion . 



Note that the decomposition in (3.2) is analogous to 



the one chosen in for the analysis of the K^^ decay. 
With this choice, the effective form factor depends only 
on the single variable t. It is convenient [|| to write it as 
the sum of two terms. 



f+{t) = .f+{t) + u\ 



(3.4) 



The first one contains the pure QCD contribution (in 
principle to any order in the chiral expansion) plus the 
electromagnetic contributions originating from the non- 
derivative Lagrangian 

£,.^0 = e^F^Z{QTQT)- (3-5) 
The explicit form of this part is given by the formula 

f+{t) = 1 + 2H^+^o{t) + HK+Ko{t) + ..., (3.6) 

where the ellipses indicate contributions of higher order 
in the chiral expansion. The meson loop function HpQ{t) 
is displayed in Appe ndix 
The second term in (3.4) represents the local effects of 
virtual photon exchange of order e^p^ , 



Ana 
327r2 V 



3 -flog 



3 log ■ 



M2 



(3.7) 



The symbol Kl2{ij) denotes the renormalized (scale de- 
pendent) part of the coupling constant K12 introduced by 
Urech in the effective Lagrangian £e2p2 describing the 
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interaction of dynamical photons with hadronic degrees of 
freedom , 15 . The couphng constants Xi , enter the 
game once also virtual leptons are taken into account . 
The coupling constant Xq (/i) is obtained from Xq (/i) after 
the subtraction of the short distance contribution B, 



^6 



(3.8) 



where 



SD 



e2 , Mj 



l-^Ew(Mp,Mz), (3.9) 



which defines also the short distance enhancement 
factor Sew{Mp, Mz) to leading order. 

In order to arrive at an infrared-finite (observable) re- 
sult, also the emission of a real photon has to be taken into 
account. The radiative amplitude M'' can be expanded in 
powers of the photon energy E-y, 



where 



-■^(0) 



and -^(0) 



(3.10) 

(3.11) 
to the non- 



Gauge invariance relates -^(.i) 
radiative amplitude M., and thus to the full form factor 
F+{t,u). Upon taking the square modulus and summing 
over spins, the radiative amplitude g ener ates a correction 



p'^iy, z) to the Dalitz plot density of (2/7). The observable 
distribution is now the sum 



p{y,z) = p'- >{y,z) +p''{y,z). 



(3.12) 



Both terms on the right hand side of this equation depend 
on the full form factor F^ and contain infrared divergences 
(from virtual or real soft photons). Upon combining them, 
the observable density can be written in terms of a new 
kinematical densi ty A i [||, and the effective form factors 
f+{t) defined in (ph. 



p{y,z)=USEw{Mp,Mz)Ai \f+{t)\' 



(3.13) 



where we have pulled out the short-distance enhancement 
factor. The kinematical density Ai is given by [|| 

A,{y,z)^Af\y,z) [l + A'^{y,z)] +AY'{y,z). (3.14) 

The function A^^{y, z) arises by combining the contri- 
butions from \Mj__^-^\'^ and r{u,ml, M^+; M^). Although 
the individual contributions contain infrared divergences, 
the sum is finite. The factor A^^{y, z) originat es fro m av- 
eraging the remaining terms of jA^'^'p [see (3.10)] and 
are infrared-finite. Note that both A^^{y, z) and Ai{y, z) 
are sensitive to the treatment of real photon emission in 
the experiment. Details on these corrections are given in 
Sect, g 

Let us finally note that, in principle, the radiative am- 
plitude generates new terms in the density, proportional to 
derivatives of form factors. These terms would only arise 
at order e'^p'^ and higher in chiral pertu rbatio n theory, and 
therefore we have suppressed them in ( ^3.13| ). 



4 Real photon radiation 

We present here in detail a possible treatment of the con- 
tribution of the real photon emission process 

TT+ip+) ^ ^"(po) e+(pe) l^eip.) liPj), (4.1) 

in complete analogy with the procedure proposed in [ pT| 
and 1^ for the analysis of the K^^ decay. To this end we 
define the kinematical variable 



{Pu + p-y)^ = (p+ ~Po~ Pef 



(4.2) 



For the analysis of the experimental data, we suggest to 
accept all pion and charged lep ton energi es wi thin the 
whole TTes Dalitz plot V given by (2.11) and (2.13), respec- 
tively, and all kinematically allowed values of the Lorentz 
invariant x defined in (^^). (The variable x determines 
the angle between the momentum of the neutral pion and 
the positron momentum for given energies Et^o , E^+ .) This 
translates into the distribution 



2i27r5 



dx 



1 f (fpi, (Pp^ 



27: J pI p% 



+ -PO-Pe 



p,-p,)^l^^^l^ 

pol 

(4.3) 



with 



a^max = + +r, 

1 



yz-f V(y2 -4re)(z2 -4r^) }. (4.4) 



In (4.3) we have extended the integration over the whole 
range of the invariant mass of t he unobserved Vi 7 system. 
The integrals occurring in (43) have the general form 



l{Pl,P2;P,M'y) 

d^q d^k 



d^^'>{P-q-k) 



1 

2^ J ko (pi • fc-^Af2/2)™(p2 • fc + Af2/2)' 



(4.5) 



The results for these integrals in the limit = c an b e 
found in the Appendix of |l^. Us ing the definition (4.5), 
the radiative decay distribution (O) can be written as 



p<'°\y,z)Ioiy,z;My) + 



327r3 



J dx^Cjn,nIni,n{Pe,P+;P+ - Po -Pe,0) 



(4.6) 
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where the infrared divergences are now confined tc 
h{y,z;M^) = 



dx 



2p+ ■ PeIl,liPe,-P+;P+ - PQ -Pe\M^) 



Mi 



Ml+ Ioa{Pe, -P+;P+ -pa-Pe;M^ 
Tnll2fi{pe,-P+]P+ -pa-Pe]M^) 



(4.7) 



The explicit form of the function /q can be found in Eq. 
(27) of ll^. (Of course, the appropriate substitutions K 
7r+ and tt — > 7r° have to be performed.) The coefficients 
Cm,n were given in Eq. (19) of Note however the mis- 
print for the values of c_i,o and ci,_2 (see Erratum of 
]l7t). As we are neglecting the contribution of the form 
factor /_. it is sufficient to consider these coefficients for 
^ = 0. 



The function introduced in (3.14) can now be re- 
lated to /q by 



loiy, z- M^) + ]-r{u, ml, M^+;M^) 



(4. 



An analytic expression of the integral occurring in (4.6) 
was given in Appendix B of [ p8[ in terms of the quantities 
U,: 

7 

^ ^m^nlm^n — ^ ^ Ui. (4.9) 





^max 



As already noticed in g, the quantity Jg{i) given in Eq. 
(A9) of |18| (which is needed for the evaluation of J/y) 
contains two mistakes: the plus-sign in the last line of 
(A9) should be replaced by a minus-sign, and 1/3™'*'^ | at 
the end of the first line of (A9) sho uld sim ply read /3™'*'^. 

The function Z\^^ introduced in (3.14) can now be ob- 
tained as 

7 



2a 



7rM\ ^ '\i=o 

7r+ i=0 ^ 



(4.10) 



5 Numerical analysis 

In the kinematically relevant range, 

ml<t< (M^+ - M^of, 



(5.1) 



the i-dependence of the effective form factor can be ap- 
proximated by the linear expansion 



/+(i) = /+(0) (1 + ^A+ 



(5.2) 



with an excellent degree of accuracy. The observable decay 
rate 

r{Tre3{j)) r{Tre3) + riTTea-,) (5.3) 

^ The right-hand side of the corresponding expression (6.7) 
in should be multiplied by 1/4. 



can now be written as 

r(7r,3(^)) = ^^ s^wiMp, Mz) \f+m^ ^(A+), (5.4) 

where 



/(A+)= dydzAiiy,z){l + —^X+ 



t 



AP 



■D 



(5.5) 



In order to extract \Vud\ we have to provide a theoretical 
estimate of the form factor at t = and the phase-space 
integral. 



5.1 Numerical estimate of /+(□) 

In the isospin limit, /+(0) coincides with the vector form 
factor at zero momentum transfer, and is thus equal to 1, 
due to the conservation of the charged isospin current. In 
the real world, all deviations from this value are therefore 
isospin suppressed. 

At one-loop acc uracy, the quantity /+ (0) obtained from 
the formula in ( 3^ ) is unambiguously determined in terms 
of the masses of the pseudoscalar mesons. It deviates from 
1 by a tiny term quadratic in the isospin-breaking param- 



eters m„ — nid and e lit 



/+(0) = 1 - 7 X 10- 



(5.6) 



Further higher orders contributions to /_|_(0) are negligibly 
small because of the aforementioned isospin suppression. 

Therefore, the theoretical prediction of /+(0) requires 
a reliable estimate of the purely electromagnetic contribu- 
tion /4_|e2p2. A numerical value for the coupling constant 
Kt^(ii) entering in (3.7) has been given by Moussallam 



Kl^iMp) = (-4.0 ±0.5) X 10 



-3 



(5.7) 



For the (unknown) "leptonic" constants we resort to the 
usual bounds suggested by dimensional analysis. 



\X^{Mp)\ <6.3x 10- 



(5.8) 



Using these numerical values, (3.7) implies 

f+\,2p2 = (4.6 ± 0.1 ± 0.4 ± 0.3) X 10-3 

= (4.6±0.5) X 10-^ (5.9) 

The three errors in the first line correspond to the un- 
certainties of Kl2{Mp), Xi and X^{Mp), respectively. For 
the final value, they have been added in quadrature. 

Despite the poor present knowledge of the "leptonic" 
constants, the uncertainty in the electromagnetic sector 
affects the final result for the effective form factor at zero 
momentum transfer by only ±0.05%, 



/+(0) = 1.0046 ±0.0005. 



(5.10) 
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Table 1. Coefficients entering the phase space integral with an associated uncertainty 





0.0 


ai 


a2 




7.375 X lO"** 


4.999 X 10"" 


1.227 X 10"^* 


a / 


7.383 X 10"* 


5.011 X 10"" 


1.230 X 10"^** 



5.2 The phase space factor 

The theoretical prediction for the slope parameter 

df+it) 



dt 



(5.11) 



t=o 



is determined by the size of the low-energy constant Lg. 
With 

Ll{Mp) = (6.9 ± 0.7) X 10"^ (5.12) 



we find 



A, 



0.037 ±0.003 



(5.13) 



The numerical coe fficients ao,i,2 entering, in the phase 
space expression (|5.5|) are shown in Table ^ Because of 
the smallness of the coefficients 01^2, the final value of 
the Dalitz plot integral /(A+) is practically insensitive to 
the exact size of the slope parameter and simply given by 
the parameter qq. The inclusion of radiative corrections 
as described in Sec. increases the value of /(A+) by only 
0.1%. 

We have also evaluated numerically oq in the case cor- 
responding to the fully inclusive one-photon decay (in- 
cluding the whole 4-particle phase space), finding no ap- 
preciable difference from the result in Table |^. 



"2 BR 



/+(0) 



(5.18) 



The present experimental precision for the branching 
ratio of the pionic beta decay cannot compete yet with the 
very small theoreti cal u ncertainty in the determination of 
Vud generated by ( 5.10| ). Using the latest value given by 
the Particle Data Group (PDG 2002) Hi, 



BR= (1.025 ±0.034) x 10-^ 

together with ( ^.10[ ), we find 

\V^^\ = 0.9675 ± 0.0160(exp.) ± 0.0005(theor. 
= 0.9675 ±0.0161. 



(5.19) 



(5.20) 



However, a substantial improvement of the experimen- 
tal accuracy is to be expected in the near future. The PI- 
BETA experiment ^ aims at measuring the branching 
ratio with a precision of about 0.5% in its current phase. 
Inserting the present preliminary result obtained by the 
PIBETA Collaboration g, 

BR== (1.044±0.007(stat.)±0.015(syst.))x 10"^ (5.21) 

we find 

\V^a\ = 0.9765 ± 0.0080(exp.) ± 0.0005(theor.) 

= 0.9765 ±0.0080, (5.22) 

to be compared with the current PDG value Q|, 

\Vud\ = 0.9735 ± 0.0008. (5.23) 



5.3 Determination of \Vud\ 



6 Conclusions 



The Kobayashi-Maskawa matrix element \Vud \ can be ex- 
tracted from the tTes decay parameters by 



\Vud\ 



7r3/2 r(^e3(.))^/' 



Gf m'(' 5ew(Mp,Mz)i/2 |/+(0)| /(A+)i/2- 

(5.14) 

We recall at this point that, according to our conve ntion 
0, the Fermi couphng constant Gp appearing in ( 5.14 ) 
has to be identified with the muon decay constant. For the 
short distance enhancement factor we use the value [|l6| 



Sew{Mp,Mz) = 1.0232, 



(5.15) 



where leading logarithmic and QCD corrections have been 
included. With the present tt^ mean life time 0, 

r^± = (2.6033 ± 0.0005) x 10"*s, (5.16) 

we finally obtain the relation 



\Vud\ - 9600.8 ^BR(7r,3(^))/ |/+(0)|, (5.17) 



The present work was devoted to the study of the pionic 
beta decay at the one loop level, with the order a radia- 
tive corrections to the amplitude included. We have been 
working within the framework of the effective low-energy 
theory of the standard model. Our ma in results in this 
respect are given by Eqs. (3^) and (3.7). 

We have discussed the influence of the various correc- 
tions on the determination of \Vud\ from a high-precision 
measurement of the pionic beta decay rate. As far as 
strong interaction corrections are concerned, the situation 
is most advantageous, since the Ademollo-Gatto theorem 
requires the deviation of /+(0) from its value 1 in the 
isospin symmetric limit to be quadratic in the meson mass 
differences - M^o and - M|.o . This results in 

a very tiny correction at one loop, ^ —7 x 10"^, and leads 
one to the expectation that higher order strong interaction 
corrections will not disturb this nice picture n. 



The situation here is very different from the case of the 
Ki3 modes, where at two loop one encounters 0[{Mfi; — M^)'^] 
counterterm contributions, which can have an influence on the 
determination of |Kis|, see e.g. the discussion in 
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Electromagnetic corrections induced by the exchange 
of virtual photons involve several unknown counterterms. 
However, naive dimensional analysis indicates that their 
contribution also remains small. For instance, they af- 
fect the extraction of |T4d| at the 0.05% level only. How- 
ever, they represent the main source of theoretical error 
at present. Putting together the short-distance corrections 
(Sew) and the long-distance corrections (to form factor 
and phase space) , we estimate an overall radiative correc- 
tion to the partial width of (+3.34 ± 0.10)% which is very 
close to other estimates [^,^. We stress that our number 
has been obtained within a completely model-independent 
framework for the long-distance corrections. 

Thus, the pionic beta decay is very close to a theorist's 
paradise, and a very precise prediction for its branching 
ratio can be obtained. It remains to be seen whether the 
experimental progresses will eventually be able to reach a 
comparable precision, and thus provide a very clean and 
accurate determination of the Cabibbo angle. 
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Appendix 

A Photon Loop Functions 

The photon loop contributions to the tTes form factor de- 
pends on the electron mass mg, the charged pion mass 
and the relativistic invariant u = (p+ — Pe)^- In or- 
der to express the loop functions in a compact way, it is 
useful to define the quantity 



X 



y - ViP 



(A.l) 



where y has been defined in (2.4). In terms of r^, X,y and 
the dilogarithm 



Li2{x) 



r^log(l-xt), (A.2) 

JO 



the functions contributing to r{u, ml, ; M^) are given 



by i 



2 log — — — ( 1 



C{u,ml,M^) = 



M2 



1 



X 



1 log2 X -f 21ogX log(l -X^)-^ 

2 b 



■log re 



+ Li2{X^) + Li2il 



X 



Li2il-X^) 



(A.4) 



ri(u,m2,M2+) = - Inre + {A - 3y)nu,ml M^^) 
r2{u, ml Ml^) ^ i (l - [ - T{u, mlM^){l - r,) 



and 



In Tr- 



T{u,mlMl^) 



-(3-y)T{u,mlM^^), (A.5) 



X 



^i-x^ 



InX 



(A.6) 



B Meson Loop Functions 

The loop function HpQ{t) 0,^2) is given by 



HpQit) 



^0 



(B.l) 



where 



(t - iSpQ) 



18(47r)2 
1 f 2SpQ - t 



12 \ ApQ 

^2[Ap{^l) + AQ{^i)] 



[Ap{n) - Aq{p)] 



with 



\{x, y, z) = + y'^ + — 2{xy + xz + yz), 



SPQ = Ml + M| , ApQ = Ml - M|, 



Ml 



Ml 

(4^)2" 



and 



JpQ{t) = 



327r2 



2 + dm. log - ^ log 

^ Ml ApQ ^ Ml 



Ml E 



PQ 



(B.2) 

(B.3) 
(B.4) 

(B.5) 



Ml 



\y\t,Ml,Ml) 



X log 



■[t + Ai/2(t,M2,M2)]2_z\ 



^PQ 



[t^X^/^t,Ml,MlW-A 



(B.6) 



The quantity Hpo(O) appearing in the evaluation of /+(0) 



is given by [11 



Ml 



hoix) = 1 + 



2x 
1 — a;' 



■ logx. 



(B.7) 
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The chiral one-loop corrections comply with the AdemoUo- 
Gatto theorem through the property 



1 



/\2 



1927r2 Fo'SpQ 



for ApQ ^ SpQ. For the theoretical determination of the 
slope parameter we need the derivative of the function 
HpQ{t) at t = given by 



dHpQ{t) 



dt 



t=o 



hi{x) 



MpMq 




2{x-l) 



(B.9) 
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